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Abstract 

(N 

^ ' It is argued that the definition of the twist-3 transverse momentum dependent correlation func- 

CN i tions must be modified if they contain Hght-Hke Wilson lines. In the framework of a simple 

^N ' spectator model of the nucleon we show explicitly the presence of a light-cone divergence for a 
specific twist-3 time-reversal odd parton density. This divergence emerges for all eight twist-3 

f^ I T-odd correlators and appears also in the case of a quark target in perturbative QCD. The 

\^ , divergence can be removed by using non-light-like Wilson lines. Based on our results we ar- 

^^ ' gue that currently there exists no established factorization formula for transverse momentum 

fH , dependent twist-3 observables in semi-inclusive DIS and related processes. 

I , Partonic correlation functions (distributions and fragmentation functions) that depend on both 

fl_j ' the longitudinal and transverse parton momentum (p^-dependent) have received considerable at- 

rS^ • tention over the past years. The most prominent examples are the (naive) time-reversal odd (T-odd) 

Sivers parton density [1, 2] and the Collins fragmentation function [3], where the status of the mere 

K^ I existence of T-odd parton distributions was clarified only relatively recently [4, 5]. A great deal of 

^ • theoretical and experimental effort has been devoted to determining these functions by characteriz- 



ing and measuring single spin asymmetries in semi-inclusive deep-inelastic scattering (DIS) [6, 7, 8] 
and azimuthal asymmetries in e"''e~-annihilation [9] (see also [10, 11]). 

The first experimental studies on single spin asymmetries in semi-inclusive DIS dealt with the 
observables Aul (unpolarized lepton beam and longitudinally polarized target) [12, 13, 14, 15] 
and Alu (longitudinally polarized lepton beam and unpolarized target) [16]. By contrast they 
are suppressed like 1/Q (with Q denoting the virtuality of the exchanged photon) with respect to 
the leading twist asymmetry Ajjt (transversely polarized target). Consequently, a parton model 
description of these asymmetries invokes (pT-dependent) twist-3 correlation functions [17, 18, 19] 
(see, e.g., Refs. [20, 21, 22, 23, 24, 25, 26, 27, 28] for related theoretical work on twist-3 single spin 
asymmetries). A complete list of p^-dependent twist-3 correlation functions was obtained only 
recently [29, 30, 31, 32, 33]. 

It is non-trivial to establish QCD-factorization for transverse momentum dependent hard scat- 
tering processes in terms of pT-dependent partonic correlators. While the twist-2 case was studied 
in some detail in the literature (see, e.g., Refs. [34, 35, 36, 37, 38]), for twist-3 observables only 
a tree- level formalism is available up to now [17, 18, 19, 39]. To some extent the difficulty is in- 
timately related to the task of finding a suitable Wilson line that ensures color gauge invariance 
of the pT-dependent correlation functions (see in particular the discussion in [40] ) . In this context 



it was emphasized already in Ref. [41], that hght-hke Wilson lines typically lead to divergences. 
At low-order in Oqcd, in calculations of twist-2 pj'-dependent correlators with light-like Wilson 
lines indeed certain diagrams show such a light-cone divergence. On the other hand, one-loop 
computations of the T-odd Sivers function defined with a light-like Wilson line yield well-behaved 
results [5, 42, 43, 44, 45, 46]. 

The main intention of this note is to delineate the presence of a light-cone divergence explicit 
in the case of twist-3 T-odd correlators with light-like Wilson lines. To this end we perform a 
one- loop calculation of the particular twist-3 T-odd parton distribution g-^ [32, 33] which appears 
in the description of the beam spin asymmetry A^u [32]. Our study is based on a simple diquark 
spectator model of the nucleon. The light-cone divergence which appears for g^ also shows up for 
all twist-3 T-odd parton distributions. Moreover, one arrives at the same conclusion when using 
a quark target in perturbative QCD. It is therefore evident that twist-3 pT-dependent correlators 
containing light-like Wilson lines are undefined. The presence of the light-cone divergence also 
implies that the frequently used tree-level formalism of Refs. [17, 18, 19, 39] for twist-3 observables 
in semi-inclusive DIS and related processes at least needs a modification. Moreover, the so-called 
subtraction formalism [48, 36, 37, 38] used for an all-order description of pT-dependent twist-2 pro- 
cesses cannot directly be applied in the twist-3 case. As a consequence, currently there exists no 
established factorization formula for processes like the longitudinal twist-3 single spin asymmetries 
Alu and Aul- 

Now, we discuss elements of the calculation of g^. The starting point is the quark-quark 
correlation function 

<^\fix,pT;P,S) = I ^^^^e'P<{P,S\^,iO)Vexp^-tg J^^ dXiv.A)iXv) 

X Vexp Ug p dX (v ■ A){Xv + o| ^»(0I^, S) 



, (1) 



where P, S, and p respectively denote the nucleon momentum, nucleon spin, and quark momentum. 
We choose a frame in which the light-cone momentum P~^ is large, and the plus-component of the 
quark momentum is given by p^ = xP^ . The two Wilson lines that appear in the definition (1) are 
determined by the vector v = {v~^ ^ v^ , Ot), with the proper light-like limit given by (0, 1, Ot)- Note 
that the twist-3 formalism of Refs. [17, 18, 19, 39] is working with light-like lines. The Wilson lines 
are of particular importance to enforce color gauge invariance in the quark correlation functions 
as well as to provide a mechanism to generate final state interactions and the associated complex 
phases necessary to characterize the T-odd parton densities [4, 5]. 

We perform a one-loop calculation in Feynman gauge for which a transverse gauge link at light- 
cone infinity [42, 47] is not relevant. The function g is defined via the correlator in (1) by means 
of [32, 33] 

2e^PTj g^ix,p^) = -P+ TT(<^ix,pT) tSs) , (2) 



where the conventions e*^^^^ = 1 and Sj, = e ^^-^ are used. Note that in (2) we take an average over 
the nucleon polarization, as the function g^ is defined for an unpolarized target. 

In the calculation the vector v is taken to be slightly off the light-cone {v^ ^ , |w+|/t;^ ^ 1) 
and space- like (f^ < 0). For our purpose here one equally well could use a time-like Wilson line. 
The light-cone divergence of g^ will show up in the light-like limit \v^\/v~ -^ 0. Non-light-like 
Wilson lines as regulators were already used much earlier in connection with a proper definition of 
pr-dependent correlation functions [41]. 
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Figure 1: One-loop diagram relevant to the calculation of g-^ for the quark in the scalar diquark 
model of Ref. [4]. [h.c] stands for the Hermitian Conjugate diagram. The thin dashed line depicts 
the cut of the correlator. In this diagram the diquark (dashed line) interacts through the exchange 
of a photon with the eikonalized quark (double line). 



We compute the correlator in the simple diquark spectator model of the nucleon used in Ref. [4]. 
In this approach the interaction between the nucleon, the quark, and the diquark is given by a point- 
like scalar vertex with the coupling constant A. We use the model with an Abelian gauge boson. 
In order to determine the distribution g to lowest non-trivial order the diagram in Fig. 1 has 
to be evaluated, where only the imaginary part of the loop amplitude (with loop-momentum /) is 
relevant. The contribution of this diagram to the correlator is given by 
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with rfi^ = x(l 



M2 + 



X "•" 1- 



and quark momentum in light-cone coordinates given 



by p 



xP^ 



p^+mi- 



-M'^{l-x) 



2P+{l-x) 

straightforward and one arrives at the expression 



in J 

e^PTj g 



Pt ) ■ Using the definition in Eq. (2) the calculation of g^ is 
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^^'^^^~"2(27r)3(p^+m2) A. 
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(4) 



The difficulty with the use of light-like lines can be understood by considering the dependence on 
the loop integration variable I in Eq. (4). If the vector v is taken to be light-like the /"-dependence 
in the propagator that is due to the Wilson line drops out. As a consequence the integral diverges 
(logarithmically) for /^ ^ and /~ -^ ±00. The crucial point is the presence of the factor l~ in 
the numerator of (4). This is in contrast to the corresponding calculations of the twist-2 T-odd 
Sivers function f^ and Boer-Mulders function hi, for which the box graph in Fig. 1 provides a 
well-defined finite result when using light-like Wilson lines. We also note that, as far as a light-cone 



divergence is concerned, the contribution of the box graph to the six twist-2 T-even p^-dependent 
parton densities is well behaved in the light-like limit, because there appears no factor of l~ in the 
numerator. Here we do not address possible complications (for T-even functions) in the case of a 
massless gluon. 

Performing the loop-integral in (4) yields the result 

I _^2\ eqCsA^ 1 (1 — x ( |f+|m^ \ 1 — 2: / 2 /x 

q {x,Pt) = -r, — TTT^T TTzS In ,^ , . „ m 2x 1 — x m — 

^ ^ '^^' 4(27r)4 (p^ + m2) \ x \v-{P+yJ a; V 2-x V2 



xpj^ 



m^ 



where we made an expansion in the (small) parameter |w^|/w^ in order to simplify the algebra. The 
first term in the brackets in (5) makes the logarithmic light-cone divergence explicit. Note also that 
the result is boost-invariant and that, due to the expansion in jf"*"!/!;", the result cannot be used 
in the limit x ^ 0. Eventually, we mention that in the case of the corresponding treatment of the 
Sivers function with a non-light-like line we recover the well-known (finite) result [5, 42, 44, 45, 46] 
in the light-like limit by putting |f+|/t;~ ^ at the end of the calculation. 

Let us also point out that the light-cone divergence can be avoided by means of a phenomeno- 
logical approach in the sense that we replace the proton-quark-diquark coupling according to 

1 - TT - «e 



i.e., we introduce a monopole form factor by hand. In Eq. (6), p denotes the quark momentum 
while A is an arbitrary mass scale. Then the expression in (4) is replaced by 



ij ±. ^2^ _ e^e^A^ A^(l - x) ^^ [ dH 

(7) 
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In (7) we use the light-like Wilson line because using this form factor results in an additional factor 
of l^ in the denominator rendering the integral finite in the potentially dangerous region /'*' ^ 
and /" -^ ±00. Although we get a finite result in this way we want to emphasize again that such 
an approach is purely phenomenological and cannot be derived using a Lagrangian of a microscopic 
theory/model. 

We end this short note by summarizing the main points and providing some additional discus- 
sion. 

• We have shown that twist-3 T-odd pT-dependent parton densities cannot be defined by using 
light-like Wilson lines. In order to make this point explicit we have presented a one-loop 
calculation of the parton distribution g , which appears for instance in the beam spin asym- 
metry Alu in semi-inclusive DIS. We arrive at the same conclusion for all eight twist-3 T-odd 
parton densities of a spin- 1/2 particle for which we pushed the calculation up to the level of 
Eq. (4). 



• 



While our calculation is based on a simple diquark spectator model of the nucleon, we find the 
same light-cone divergence also in the case of a quark target in perturbative QCD. Therefore, 
the main conclusion we have drawn is not based on some artefact of a simple model. 

With regard to possible light-cone divergences from Fig. 1, all eight twist-2 p^-dependent 
parton densities are well-behaved in the light-like limit. 

• Although we expect that the same divergence problem arises from the graph in Fig. 1 for T- 
even twist-3 parton densities, at present we cannot exclude that there appears a cancellation 
of divergent terms once all one- loop diagrams are taken into account. 

• For the T-odd functions we considered here the light-cone divergence can be regularized 
by using non- light-like Wilson lines (see, e.g., also Ref. [41]). We introduced such lines by 
hand. It would be very interesting to investigate (as a first step) if a tree-level factorization 
formalism in the spirit of Refs. [17, 18, 19, 39] can be established which automatically provides 
definitions of twist-3 py-dependent correlators with non-light-like Wilson lines. 

• In the literature a subtraction method has been used (see, e.g., Refs. [48, 36, 37, 38]) in order 
to formulate all-order factorization for transverse momentum dependent twist-2 observables. 
In this formalism also correlation functions with non-light-like lines appear. However, in 
particular when performing one-loop calculations of an observable, the dependence on various 
non-light-like directions drops out by construction, as it should be. In contrast, using this 
formalism in the case of, e.g., the longitudinal asymmetry A-'^jj for jet production in DIS, 
which exclusively is given by g^ [32], leads to a result depending on the a priori arbitrary 
ratio \v^\/v^. It is quite possible that the subtraction formalism can be generalized such 
that it can cope with transverse momentum dependent twist-3 observables. But this requires 
further studies. 

• In a very recent manuscript [49] also a result for g in the scalar diquark model has been 
provided. Rather than starting from an operator definition, this result is extracted from 
the calculation of the beam spin asymmetry Am in semi- inclusive DIS. By doing so and 
assuming factorization in the spirit of Refs. [17, 18, 19, 39] one arrives at an expression for 
g^ which does not show the divergence discussed in our note, because the result for the full 
asymmetry Am is finite in the diquark spectator model [30, 31, 49]. However, apparently 
the factorization assumption is not valid for this twist-3 process, and g^ obtained in Ref. [49] 
is not the corresponding twist-3 distribution function defined through Eqs. (1) and (2). 

Acknowledgements: We thank Gary R. Goldstein for useful discussions. The work of M.S. has 
been supported by the Graduiertenkolleg "Physik der Elementarteilchen an Beschleunigern und 
im Universum." The work has also been partially supported by the KOSEF (Korea Science and 
Engeneering Foundation), by the Verbundforschung (BMBF), and by the Transregio/SFB Bochum- 
Bonn-Giessen. This research is part of the EU Integrated Infrastructure Initiative Hadronphysics 
Project under contract number RII3-CT-2004-506078. 

References 

[1] D.W. Sivers, Phys. Rev. D 41, 83 (1990). 
[2] D.W. Sivers, Phys. Rev. D 43, 261 (1991). 
[3] J.C. Collins, Nucl. Phys. B396, 161 (1993). 



[4] S.J. Brodsky, D.S. Hwang, and I. Schmidt, Phys. Lett. B 530, 99 (2002). 

[5] J.C. Collins, Phys. Lett. B 536, 43 (2002). 

[6] A. Airapetian et al. [HERMES Collaboration], Phys. Rev. Lett. 94, 012002 (2005). 

[7] V.Y. Alexakhin et al. [COMPASS Cohaboration], Phys. Rev. Lett. 94, 202002 (2005). 

[8] M. Diefenthaler, arXiv:hep-ex/0507013. 

[9] R. Seidl et al. [Belle Collaboration], Phys. Rev. Lett. 96, 232002 (2006). 

[10] W. Vogelsang and F. Yuan, Phys. Rev. D 72, 054028 (2005). 

[11] A.V. Efremov, K. Goeke, and P. Schweitzer, Phys. Rev. D 73, 094025 (2006). 

[12] A. Airapetian et al. [HERMES Collaboration], Phys. Rev. Lett. 84, 4047 (2000). 

[13] A. Airapetian et al. [HERMES Collaboration], Phys. Rev. D 64, 097101 (2001). 

[14] A. Airapetian et al. [HERMES Collaboration], Phys. Lett. B 562, 182 (2003). 

[15] A. Airapetian et al. [HERMES Collaboration], Phys. Lett. B 622, 14 (2005). 

[16] H. Avakian et al. [CLAS Collaboration], Phys. Rev. D 69, 112004 (2004). 

[17] J. Levelt and P.J. Mulders, Phys. Lett. B 338, 357 (1994). 

[18] P.J. Mulders and R.D. Tangerman, Nucl. Phys. B461, 197 (1996) [Erratum-ibid. B484, 538 
(1997)]. 

[19] D. Boer and P.J. Mulders, Phys. Rev. D 57, 5780 (1998). 

[20] K.A. Oganessian, H.R. Avakian, N. Bianchi, and A.M. Kotzinian, arXiv:hep-ph/9808368. 

[21] A.V. Efremov, K. Goeke, M.V. Polyakov, and D. Urbano, Phys. Lett. B 478, 94 (2000). 

[22] E. De Sanctis, W.D. Nowak, and K.A. Oganessian, Phys. Lett. B 483, 69 (2000). 

[23] B.Q. Ma, I. Schmidt, and J.J. Yang, Phys. Rev. D 63, 037501 (2001). 

[24] A.V. Efremov, K. Goeke, and P. Schweitzer, Phys. Lett. B 522, 37 (2001) [Erratum-ibid. B 
544, 389 (2002)]. 

[25] A.V. Efremov, K. Goeke, and P. Schweitzer, Phys. Rev. D 67, 114014 (2003). 

[26] F. Yuan, Phys. Lett. B 589, 28 (2004). 

[27] L.P. Gamberg, D.S. Hwang, and K.A. Oganessyan, Phys. Lett. B 584, 276 (2004). 

[28] M. Anselmino, M. Boglione, U. D'Alesio, A. Kotzinian, F. Murgia, and A. Prokudin, Phys. 
Rev. D 71, 074006 (2005). 

[29] K. Goeke, A. Metz, P.V. Pobylitsa, and M.V. Polyakov, Phys. Lett. B 567, 27 (2003). 

[30] A. Afanasev and C.E. Carlson, arXiv:hep-ph/0308163. 

[31] A. Metz and M. Schlegel, Eur. Phys. J. A22, 489 (2004). 

6 



[32] A. Bacchetta, P.J. Mulders, and F. Pijlman, Phys. Lett. B 595, 309 (2004). 

[33] K. Goeke, A. Metz, and M. Schlegel, Phys. Lett. B 618, 90 (2005). 

[34] J.C. Collins and D.E. Soper, Nucl. Phys. B193, 381 (1981) [Erratum-ibid. B213, 545 (1983)]. 

[35] J.C. Collins, D.E. Soper, and G. Sterman, Nucl. Phys. B250, 199 (1985). 

[36] X. Ji, J.P. Ma, and F. Yuan, Phys. Rev. D 71, 034005 (2005). 

[37] X. Ji, J.P. Ma, and F. Yuan, Phys. Lett. B 597, 299 (2004). 

[38] J.C. Collins and A. Metz, Phys. Rev. Lett. 93, 252001 (2004). 

[39] D. Boer, P.J. Mulders, and F. Pijlman, Nucl. Phys. B667, 201 (2003). 

[40] J.C. Collins, Acta Phys. Polon. B 34, 3103 (2003). 

[41] J.C. Collins and D.E. Soper, Nucl. Phys. B194, 445 (1982). 

[42] X. Ji and F. Yuan, Phys. Lett. B 543, 66 (2002). 

[43] D. Boer, S.J. Brodsky, D.S. Hwang, Phys. Rev. D 67, 054003 (2003). 

[44] L.P. Gamberg, G.R. Goldstein, and K.A. Oganessyan, Phys. Rev. D 67, 071504 (2003). 

[45] A. Bacchetta, A. Schaefer, and J.J. Yang, Phys. Lett. B 578, 109 (2004). 

[46] K. Goeke, S. Meissner, A. Metz, and M. Schlegel, Phys. Lett. B 637, 241 (2006). 

[47] A.V. Belitsky, X. Ji, and F. Yuan, Nucl. Phys. B656, 165 (2003). 

[48] J.C. Collins and F. Hautmann, Phys. Lett. B 472, 129 (2000). 

[49] A.V. Afanasev and C.E. Carlson, arXiv:hep-ph/0603269. 



